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Abstract 

Exact expressions for probability densities 
of conjugate pair separation in euclidean 
isometries are obtained, for the cosmic crys- 
tallography. These are the theoretical counter- 
parts of the mean histograms arising from 
computer simulation of the isometries. 
For completeness, also the isometries with 
fixed points are examined, as well as the 
orientation reversing isometries. 



1 Introduction 

Various methods have been proposed to in- 
vestigate the shape of the universe, and cosmic 
crystallography (CC) is one of them [|T|. As 
CC ponders, if the universe is multiply con- 
nected (MC) then multiple images of a same 
cosmic object may be seen in the sky. The 
separations between these images are corre- 
lated by the geometry and the topology of the 
spacetime; so if one selects a catalog of n ob- 
served images of various cosmic objects and 
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performs a histogram of the n(n — l)/2 sepa- 
rations between these n images, then the ex- 
isting correlations must somehow show up. 

It was recently examined @] in what re- 
spects the histogram for a multiply connected 
observed universe should differ from that of a 
simply connected (SC) one, with same geom- 
etry and radius. It was found that each isom- 
etry of the MC universe individually imprints 
either a small localized deformity on the his- 
togram of the SC universe 0], or a sharp spike 
if the isometry is a Clifford translation ||. 

Since each observed universe model has a 
unique pair separation normalized his- 
togram, a strategy to unveil the cosmic topol- 
ogy seems straightforward: one should com- 
pare the histogram obtained from observa- 
tional astronomy with histograms obtained 
from computer simulated universe models 
with prescribed topologies. 

Both types of histograms (observational and 
simulated) are infected with statistical noises, 
and methods to reduce these noises are de- 
sirable. A suggestion was made, to replace 
the histogram related to the SC component of 
the simulated model by an exact continuous 
probability density function. For each geome- 



try with constant curvature the corresponding must be obeyed, 
function was then derived [[| ; however, appro- 
priate functions were still lacking, to replace 
the histograms related to each isometry com- 
ponent p. 

In the present study we derive some of these 
functions, namely the normalized pair sepa- 
ration probability densities for the euclidean 
isometries. Following the prescriptions of ref. 
0, if we now merge these new functions with 
the functions already given in M then we ob- 
tain noiseless normalized probability densities 
more suitable for comparing with the normal- 
ized observational histograms. For complete- 
ness we also examine the euclidean isometries 
with fixed point, as well as the orientation re- 
versing isometries. 



It is clear that the balls B and B g need inter- 
sect, so the separation m between their centers 
C and C g has to satisfy 



m < 2a . 



(2.2) 



The intersection B fl B g is a rotationally sym- 
metric solid lens whose diameter, thickness, 
and volume are (see figure 1) 



D 
T 



V Aa? — m 2 , 
2a — m, 

— (2a - m) 2 (4a + m) . 



(2.3) 
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2 The euclidean isometries 

The isometries of the euclidean space E 3 are 
the (pure) translations t, (pure) rotations u>, 
screw motions tu, (pure) reflections e, and 
glide reflections et. While t, u, and tu are ori- 
entation preserving, e and et reverse the orien- 
tation of E 3 ; and while t, too, and et act freely 
on E 3 , uj and e have fixed points. We study all 
five isometries for the sake of completeness, al- 
though cosmic crystallography is presently in- 
terested on the orientation preserving isome- 
tries without fixed points, namely t and tu 
only. 

Our general approach to investigate an 
isometry g in cosmic crystallography is first 
consider a solid ball B with radius a, then ap- 
ply g to B thus producing a new solid ball B g , 
next consider the set of points P G B whose 
corresponding transformed points P g are also 
in B, finally randomly select a pair (P, P g ) 
and ask for the probability V g (l)dl that their 
separation lies between the values I and I + dl. 
The normalization condition 



a 
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j Q V^(l)dl = l 



Figure 1 The solid balls B and B g both with 
radius a intersect in a solid lens with diameter 
D and thickness T. □ 

3 Translations, screw mo- 
tions, and rotations 

Translations 

If the euclidean space is subjected to a nonzero 
translation t, the probability density that a 
point P be displaced a distance I is clearly 



(2.1) 



V t {l)=6(l-t) 



(3.1) 
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the Dirac delta; the density satisfies the nor- 
malization condition fl2.ll) , and does not de- 
pend on the ball B. 

Basics on screw motions of a solid ball 

In E 3 , imagine a straight line £ placed 
at a distance b from the center C of the 
solid ball B with radius a; £ and B may 
a), be tangent (6 = a) 
{b > a); see figure 2. 
£ 

(b < a) 

£ 

(b=a) 
£ 

(b > a) 



intersect (6 < 
or be disjoint 




Figure 2 Relative positions of a solid ball B 
with radius a and a straight line £ at a dis- 
tance b from the center of B. □ 

Now consider a screw motion g of the ball, 
with nonzero translation t parallel to the line 
£ and nonzero rotation uj around the line; for 
our purposes the senses of t and u are irrele- 
vant, so for definiteness and simplicity we as- 
sume t > and < uj < ir. 

The separation m between the centers C 
and C g of the intersecting balls B and B g is 



m 



t 2 + 46 2 sin 2 cu/2, 



(3.2) 



and the condition m < 2a implies the con- 
straint 



t 2 + 4 6 2 sin 2 u;/2 < 4a 2 



between the four independent parameters a, 
b, t, and uj; see figure 3. 




sin co/2 



Figure 3 A solid ball (not displayed) with 
center C rotates uj around the line £ (the 
z axis) and translates t parallel to the axis, 
eventually reaching the position centered at 
C g ; both centers C and C g are at a distance 
b from the axis, and their separation m must 
be smaller than the sum 2a of the two radii in 
order that the balls intersect. □ 

We next introduce the auxiliary variable r 
and a function Qg{r) which help simplify our 
study. Call r the separation between the point 
P (also P g ) and the axis £ (see figure 4); it 
is related to I (the separation from P to P g ) 
through 



f = t 2 + 4r 2 sin 2 cu/2 ; 



(3.4) 



for each isometry (t, uj) this is a bijective rela- 
tion between I and r, so the probability den- 
sity Vg(l) of finding a g-pair with mutual sep- 
(3.3) aration I parallels the akin probability den- 
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sity Qg(r) of finding a g-pair whose members of the volume V B of the intersection B H B, 
(both in B) are at a distance r from the axis eq. 
of the motion: 



9 ■ 



V*{l)dl 



Q B Jr)dr. 



(3.5) 




Figure 4 Under the screw motion g with axis 
£ , translation t , and rotation u , a point P 
separated r from L moves to the position P g 
at the distance I given by eq. (|3.4|). □ 



Clearly the normalization condition 

Q B Ar)dr = 1 (3.6) 



is also to be satisfied. Once we obtain Q B {r) 



and since from 



we have 



Vi 2 ^ 

r = T i 

2 sinu;/2 

then we will finally compute 



(3.7) 



(3. 



To obtain Q B (r) we first consider a suffi- 
ciently long cylinder C r with radius r and axis 
C. A short reflection gives that the probabil- 
ity density Q B {r) is linearly proportional to 
the area S B {r) of the intersection B fl B g fl C r ; 
the coefficient of proportionality is the inverse 



Q B g {r) 



area{B fl B g fl C r 
wo/(£n£ 9 ) 



(3.9) 



To find the area S B {r) we first assume the line 
C (the axis of the screw motion) along the 
cartesian 2-axis; the cylinder C r with axis C 
and radius a then has equation 



2 i 2 

x + y 



(3.10) 



Still without loss of generality we take the cen- 
ter C of the solid ball B with radius a at the 
cartesian position (b, 0, 0); the points of B 
then satisfy 



{x - bf 



y 



2 ^ 2 

z < a . 



(3.11) 



Finally the center C g of the solid ball B g is at 
the cartesian position (b cosu, b sinu, t), so 
the points of B g satisfy 

(x — bcosu) 2 + (y — bsinu) 2 

+ (z-t) 2 <a 2 . (3.12) 

S B (r) is then the area of the surface whose 
points (x,y,z) satisfy ( |3.10|) , ( |3.11|) , and 
( |3.12|) simultaneously. To visualize this sur- 
face B fl B g fl C r we first consider the surface 
B fl C r , then the similar surface B g fl C r , and 
finally the combined intersection (i3 fl C r ) fl 

(fi 5 nc r ). 

The form of the intersection B fl C r depends 
on the relative values of a, b, and r (see figure 
5): 

(i) if < r < a — b then all generatrices of C r 
penetrate into B: the intersection is a topo- 
logical ring (an annulus, a compact cylinder, 
a disc with a hole in it); 

(ii) if | a — b\ < r < a + b then only a part 
of the generatrices of C r goes through B: the 
intersection is a topological disc; 
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(iii) if a + b < r, or if < r < b — a, then there 
is no intersection. 



(see figure 6) 




Figure 5 (A) The solid sphere B and the 
sufficiently long cylindrical surface C r , in per- 
spective; in the other drawings the line of sight 
is the vertical; 

(i) the intersection BnC r is a topological ring 
when < r < a — b (loosely saying, a is too 
large); 

(ii) B and C r intersect in a topological disc 
when \a — b\ < r < a + b (then a, b, and r may 
form a triangle); 

(iii) B and C r do not intersect when a + b < r 
nor when < r < b — a (loosely saying, either 
r or b is too large). □ 

Since the intersections B C\C r are drawn on 
the cylinder C r itself, we use the cylindrical 
coordinates (f> (the azimuthal angle) and z (al- 
titude) to visualize them. A ring-like intersec- 
tion B fl C r is then bounded by the two curves 



z = ±v a 2 — b 2 — r 2 + 2br cos i 



(3.13) 



where (f> G (— 7r, tt); the positive curve oscil- 
lates between the extreme values 



^max 



a 2 — {b — r) 



^2 _ ( 6 + r )2 _ ( 314 ) 



2 


^max 


1 






z min 


-71 o 













Figure 6 The ring-like intersection £> fl C r , 
eq. (|3.13|) , when a = 4 , b = 19/5, and r = 
1/10 . The boundaries at <fi = —it and <fi = n 
are identified. The upper and lower bound- 
aries are not sinusoidal. □ 

In the disc-like intersection BnC r (see figure 
7) the surface is again read from ( |3.13| ), but 

now (j) E [-(f> m ax,<t>max) with 



COS 



2br 



while z max is still given in ( |3.14| ), z min 



(3.15) 




now. 



z 
0.5 




I ° 
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Figure 7 The disc-like intersection B D C r , 
eq. (|3.13|) , when a = 1 and b = r = 1 / ; it is 
an oval (not an ellipse) centered at <fi = z = , 



with 



7r/2 and z„ 



1. 



□ 



A disc-like intersection B g fl C r is a g- 
transported copy of the corresponding inter- 
section BnC r , parallelly dragged on the cylin- 
der: while BC\C r is centered at (<fi — 0, z — 0), 
the center of £> 9 fl C r is at (0 = ui, z — t). 

Similar statements are true for the ring-like 
intersections B g fl C r . By inspection it is then 
easy to visualize the form of the combined in- 
tersection B R B g fl C r in both disc-like and 
ring-like cases. Since has a cyclic charac- 
ter, the nonempty intersection of two discs on 
C r can be either 1 or 2 discs, while two rings 
intersect in either 1 or 2 discs, or 1 ring (see 
figures 8 and 9). 



Figure 9 Instance of disconnected intersec- 
tion B fl B g fl C r (dashed areas) when B fl C r 
(and Bg fl C r ) are disc-like; here a = 3.5 , b = 
r = 2,t = 0.7, and cu = 14tt/15. □ 

We are now ready to evaluate the area Sg(r) 
of the combined intersection BnB g nC r ; given 
the fixed values of a, b, t, and u, then for each 
value of r we need to integrate the differential 



dSf(r) 



Az.rc 



(3.16) 



where Az(r, <fi) is the varying height of the in- 
tersection, and where the limits of integration 
in may vary with r (see eq.( [3.15|) ). 

As is evident from figures 8 and 9, the com- 
bined intersection is always bounded by one of 
the two curves 



zx{<f>) = ^a 2 -b 2 -r 2 + 2br cos (f) , (3.17) 
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Figure 8 Instance of disconnected intersec- 
tion B fl Bg fl C r (dashed areas) when B fl C r 
(and B g flC r ) are ring-like; here a = 4.505 , b = 
2 , r = 2.5 , t = 5.3 , and uj = tt . □ 




z 3 ((j))=t + ^a 2 -b 2 -r 2 + 2brcos{<f)-uj), (3.18) 
from the upper side, and one of the two curves 
z 2 {4>) = -z x {4>) , (3.19) 

Zi ( ( p)=t-^a 2 -b 2 -r 2 + 2br cos(0-cj), (3.20) 

from the lower side. Further, we only 
have nonempty combined intersection when 
zi(4>) > z&_{4>). We then obtain 



x 



rrun[Zi, z 3 ) 



Q{Z! - Z A ) 

- max(z2, za) 



(3.21) 



with 4>max given by ( P-15| ) for disc-like intersec- 
tions, and (pmax = ft for ring-like intersections; 
is the step function with values and 1. 

Elliptic integrals 

As is seen from eqs. ( |3.17|) - (p.21|) , the area 
Sg{r) involves terms of the form 



A(r) 



a(r 



a 2 — b 2 — r 2 + 2br cos 



(3.22) 
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with cos(0 — lo) sometimes replacing cos0 in two separate discs regime for l/v8 < r < 1/2 



the integrand. To evaluate these terms we de- (l < / < J3/2), and is empty for larger r and 
fine i 

J 2 ( r ) — a 2 — (b — r) 2 ^ e ^ rs ^ cons ^ er the ring-like regime: pre- 

liminarly find, from (|3.23), 



k 2 ( r ) — i f2(„\ _ 3 | /' .2 



(3.23) 

after noting that the condition f 2 (r) > is _ (3.27) 



satisfied by all nonempty intersections B fl C r . k : (r) 



72 



r 



f 



Then introduce the (incomplete) elliptic inte- 
gral of 2 nd kind and the area ^Tj) of the intersection B n £ 5 D 



c 

E( % k)= r^zE^dx, (3.24) 



VI -x 2 ~ S®(r)=r(8fE(k)-nV2) ; (3.28) 

and finally have the integral ( |3.22| ) expressed 



9 

now since 



A(r) = 2f(r)E( S in 1 -a(r),k(r)), (3.25) r = - 1/2, dr/dl = Z/(4r) , (3.29) 

where / = V? 2 " and jfe = Vl 2 " . When 7 = 1 the probability density (jp) is, for l/y/2 < 



in ( p.24p , or equivalently a = 7r in ( |3.22 



we 



have the complete elliptic integral of 2 kind , 

E(k) = E(l, k) ; (3.26) AV g V 7 

the functions E(k) only appear in the ring-like where the volume of B H £ 9 is 
intersections. 

Which of Zi or Z3 is minimum in ( |3.21| ), and = 7r (4v / 2 — 5)/12 . (3.31) 
which of Z2 or z 4 is maximum and whether 

21 > z 4 or ^ < 2 4 , usually depends on the Clearly /(r) and k(r) in O) need consider 

angle 0; clearly also the values of the param- the dependence r(l) given in flp9|). 

eters a,6,t, and w are relevant, well as the To describe the two-discs regime we further 



radius r. define 

A sample screw motion with b 7^ g(r) = y / l/(2r) 2 — 1 , 

We consider the case with a = t = 1/ a/2, 6 = 1 

1/-/8, = 7r; the intersection BnC r is then #±( r ) = \ -, (3.32) 



ring-like for < r < l/\/8 (equivalently ' 

l/\/2 < I < 1 ), is disc-like for 1/Vo < . r ,•, ^rrrnA 
' „ , r- , ' ' , . then obtain tor the area (3.21) 

r < 3/Vo (or 1 < / < yo), and is empty 

for larger values of r and I. The combined in- ^ = grf ( E( k) - E{q_ k) 

tersection B fl B g fl C r has a ring-like regime 9 ^ V ^ ' ^ o 

for < r < 1/V8 (1/a/2 < / < 1) , then a + \/2r(2cos _1 5 - vr) , (3.33) 
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and finally, for 1 < I < y3/2 , the probability 
density 



I 



4V* 



8f(E(g + ,k)-E(g_,k) 



+ V2(2 cos -1 



7T 



(3.34) 



In ( p.34| ) we clearly need replace r by its value 
r(l) in (|3.29| ) in the functions f,g,g + ,g_, and 
k. In figure 10 we reproduce the function 
Vg(l) for this sample screw motion, for the 

whole interval l/y/2 <l < Js/2. 



provided r < a; since we assumed t > 
0, we have Tnin(zi,Zs) = V 'a 2 — r 2 and 
max(z2, z±) = t — \J a 2 — r 2 , while the condi- 
tion z\ > Z4 implies 

r < r max = ^a 2 - t 2 /A . (3.37) 

The of the combined intersection is then 

Sf(r) = 2nr{2Va 2 - r 2 -t) [b = 0] (3.38) 

whenever positive, otherwise it is zero. Finally 
the probability density is, when 6 = 0, 

dr S*(r) 



9{> dl VP 



6/(2v^ 



t) 



0.7 



csc 2 o;/2 (3.39) 

if positive, otherwise it is zero; in ( |3.39| ) one 
clearly has to substitute r for its expression 
eq. (|3~7j) . As a matter of fact we have Vg(l) ^ 
1 1.2 1 when b = only for l m in < I < Imax, where 



Figure 10 The probability density for 
the screw motion of the solid ball B with a = 
2b = t= 1/V2, and uo = vr. For I G [1/V2, 1] 
we use the ring-like equation (|3.30|), while for 



Imax = \Jt 2 cos 2 uj/2 + 4a 2 sin 2 cj/2 , (3.40) 
and we have in these extreme limits 



I G [1 , y3/2] the two-discs equation (|3.34| ) is 
used. The integrated area is 1. □ 

The special cases with 6 = 

In these cases the axis C of the screw motion 
contains the centers C and C g of the solid balls 
B and B g ; the intersections B H C r and B g fl C r 
are ring-like and have constant width, as well 
as the combined intersections Br\B g r\C r . From 

gives 



V B g (l min ) = Qt{2a - t)-\4a + t^csc 2 w/2 . 



V B {] 

' g \ b max ) 







(3.41) 



a sample graph of V g (l) for screw motions 
with b = is given in figure 11. 



B 
V g (l) 



72|) we have m = t, then (P^3 

7T 
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Vf = —{2a - t) 2 {Aa + t) [b = 0] (3.35) 



provided t < 2a; this sine qua non is assumed 
whenever 6 = 0. The equations (|3.17|) -( |3~20"D 
simplify to 



1.4 1 



Z\ = —Z2 = Va 2 — r 2 , 

^3 = t + Z\ , Z4 = t — Z\ 



Figure 11 The probability density V g (l) for 
(3.36) the screw motion of the solid ball B with 6 = 
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0, a = 1/ v^2, t = 1 , and c<j = 7r , eq. ( |3.39|) . 
The integrated area is 1. Note that l min = 1 , 
l max = V2 , and Vf(l) = 6(3 + y/2)/7 ~ 3.8 . 
□ 

Rotations 

A rotation is the limit of a screw motion when 
the translation tends to zero, so it is formally 
described by simply setting t = in the ap- 
propriate preceding equations. In particular, 
eq. ( |3.4j ) becomes I = 2rsmu>/2. Rotations 
of the solid sphere B are classified into two 
categories: 

(i) if a < b then B is exempt of fixed points; 

(ii) if b < a then the axis £ traverses B, which 
has fixed points. 

In category a < b the intersection condition 
( |3.3| ) clearly must be satisfied (see figure 12), 



b sin <jjj 2 < a 



(3.42) 




All intersections BnC r (and B g nC r ) are topo- 
logical discs. Since Z\ > and = —Z3 in 
rotations, the area ( 3.21 ) has the simpler ex- 
pression 



S?(r) 



2r 



min(zi, z 3 )d(p . (3.45) 



When b smcj/2 < a < b and r m — D/2 < r < 
r m + D/2 we obtain 

Sf(r) = 8rf\E(l/k, k) - £(sin u/A, kj , (3.46) 



with f(r) and k(r) given in ( |3.23| ). The proba- 
bility density (|3.8|) is easily written if we take 



in succession I = 2rsino;/2, then Qg(r) = 
Sg(r)/V g B with (|3~46| ), finally Vf in Q with 
m = 2 b sinu;/2 ; in figure 13 a sample graph 
of Vq(l) for rotation with a < b is given. 




Figure 12 The solid lens B H B g when the 
solid ball B rotates u around the ^-axis and 
6sinu;/2 < a < b; the lens lays between the 
radial positions r min and r max given by eq. 

The lens B fl B g lays between the radial po- 
sitions 

Tmin '"m D/2 , T max T m -\- D/2 , 

r m = b cosu/2, (3.43) 

where r m is the radial coordinate of the center 
of the lens and D is its diameter (12.31) with 



m 



2 b sin oj/2 . 



(3.44) 



Figure 13 The probability density V g (l) for 
a rotation u of the solid ball B with radius 
a < b . Here a = 1 , b = 2 , and u = tt/6 . The 
integrated area is 1. □ 
In category b < a the intersections BnC r are 
either ring-like or disc-like, and the combined 
intersections B fl B g fl C r can be of three types 
according to the relative values of a, b, u, and 
r (see figure 14): 

(i) a ring, when < r < a — b ; then 

S^(r) = 8rf [2E{k) 

-E{sinuj/4, k) -E{cosu/A, k)] ; (3.47) 
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(ii) a pair of discs, when a — b < r < D/2—r m ; given, for the entire interval 0<r <D/2+r r , 
then 



0.3 



Sf(r) = 8rf[2E(l/k, k) 

-E{wa.u/1, k) -E(cosu/A, k) 1 ; (3.48) 



(hi) a disc, when D/2 — r m < r < D/2 + r m ; 
then 

Figure 15 The probability density for a 

(pure) rotation of the solid ball B with radius 
a > b. Here a = 2,6 = 1.3, and u = n/3. 
Sf(r) = 8rf[E(l/k,k)-E(smuj/4,k)]. (3.49) The integrated area is 1. □ 





two discs 



ring 

Figure 14 The solid ball B with center C and 
radius a > b rotates u around the 2-axis and 
gives the new ball B g with center C g . The 
solid lens B fl B g is intersected by the cylindri- 
cal surface C r with radius r and axis along the 
z-axis. The intersection B fl B g fl C r is a ring 
if r is small (0 < r < a — 6) , is a pair of discs 
for a — b < r < D/2 — r m , and is a single disc 
when r is larger (D/2 — r m < r < D/2 + r m ). 
□ 

In figures 15 and 16 sample graphs of the 
probability Vg(l) for rotation with b < a are 




Figure 16 The probability density V g (l) for a 
(pure) rotation of the solid ball B with radius 
a > b . Here a = 1 , b = 5/7 , and u = 2 sec _1 5. 
The region where each of the three expres- 
sions (|3.47|) - (|3.49|) for S g (r) is used is indi- 
cated. The integrated area is 1. □ 

Rotations with 6 = 

In the category b < a the special cases 6 = 
deserve a few words. The solid spheres B and 
B g coincide, and we easily find that 



I = 2r sin u/2 , S 
Vf = 4vra 3 /3 . 



47rrv / a 2 — r 2 , 
(3.50) 
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The probability density is then 



3/\/4a 2 sin 2 cj/2-Z 2 
° { ' 8 a^m 3 u/2 



(3.51) 



between the three independent parameters 
a, 6, and t; a fortiori 6 < a, so both balls in- 
tersect the plane X ; see figure 18. 



whenever < I < 2a sin uo/2 ; a sample graph 
is given in figure 17. 




Figure 17 The probability density V g (l) for 
a (pure) rotation with 6 = 0. Here a = 1 , and 
= 7r/6 . The integrated area is 1 . □ 




Figure 18 The solid ball B with center C and 
radius a is first reflected on the plane x = 
a distance b < a apart, then translated t 
upwards. □ 

Now randomly choose a point P = (x, y, z) 



of B, such that P g 
clearly P G B(~)B„ 



= (—x, y, z + t) is also in i3; 
Z being the separation from 
P to P g , we ask for the probability density 
V® (I) as described in sec. 2; we readily find 



4 Reflections and glide re- 
flections 

In E 3 , let a solid ball B with radius a have 
its center C at the cartesian position (6,0,0). 
Next consider a glide reflection g on the plane 
x = (the plane Aq) with nonzero translation 
t in the direction +z. The center C 9 of the 
glide reflected ball B g is at the cartesian posi- 
tion (— 6, 0,t), and the separation m between 



C and C g is 



m 



V46 2 + 1 2 



The condition of intersection B H B g 
plies the constraint 



46 2 + t 2 <4a 2 



(4.1) 



lm- 



(4.2) 



V B g(l) 



area(B n B g n X x 



vol(BnB g 



where is one of the two planes 



1 



x = i-VZ 2 ^ 2 ; 
2 



(4.3) 



(4.4) 



the multiplying factor 2 in (|4.3|) accounts for 
these two possibilities, and the volume V g B is 



given in fl2.3|) . 

The intersections B n £> g n are of two 
types, depending on whether or not the plane 
X x intersects the equator £ of the lens B fl B g 
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(see figure 19): 



(ii) a circular disc with radius 




c 



M 



x 



A = 



^rnax 
(i) bicircular 




a- b= x, 



max 



(ii) circular 




Figure 19 The solid lens B D B g and the two 

types of intersections B H B g H X x when the 
isometry is a glide reflection. In the equator 
£ of the lens, M is the point placed farthest 
from the reflector plane x — . In (A), M is 
below the center C g of the glide reflected ball 
B g , then x max = x M = A given in (£0]); only 
bicircular sections (i) are possible. In (B), M 
is above C g , then x max = a — b > % and two 
types of vertical cross-sections can occur: (i) 
bicircular, and (ii) circular. □ 

(i) a bicircular disc, when X X C\S ^ (ft; the disc 
is enclosed by two circles with unequal radii 



a 2 — (b =F x) 2 



(4.5) 



r±(xj = 
It has area 

oBf \ 2 -lh- i 2 -1^4- ^ //I p\ 

d" (x) =r_ cos hr, cos , (4.6) 

y r_ r + 2 

whenever |x| < A , where 

h±(x) = t/2±2bx/t, (4.7) 



t 



A 2^6 2 + t 2 /4 



- 1 



8(x) = 2 y / (l + 46 2 /t 2 )(A 2 



(4.t 



p = v /a 2 + (6+|x|) 2 , (4.9) 

which appears only when 

g = 46(a-6) -t 2 >0, A<|x| < a-b; (4.10) 

the area of the disc is clearly irp 2 . 

Collecting together preceding terms we fi- 
nally obtain 



V B Al) 



l 



l 



VP~^v g B 



0(A-\x\)S*(x) 



+Q(q)Q(\x\ - A)Q(a-b- \x\)np 2 , (4.11) 



valid for t < I < l max where 



.,„,„■- V* 2 + 4(a- ft) 2 9(g) 
at 



+ 



6 2 + t 2 /4 



(4.12) 



A graph of Vf(l) is given in figure 20. 




Figure 20 The probability density V g (l) for 
a solid ball under a glide reflection, eq. (|4.11| ). 
Here a = 2 , b = 1 , and t = 0.25 . The regions 
where the sections BV\B g V\X x are bicircular or 
circular are displayed. The function diverges 
when I — > t , nevertheless the integrated area 
from t to lmax is finite. □ 

Reflections 

(Pure) reflections are glide reflections whose 
translation is t = . 
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Bo 




b 111 

Figure 21 The solid lens B H B g when the 



isometry g is a (pure) reflection on the plane 
x — . The sections i 
with variable radius p . 



a; = . The sections B H B g H X x are circles 

□ 



As is evident from figure 21, the probability 
density Vg(l) is proportional to the of a 
disc with radius 



p=y/a?-(b + l/2)< 



(4.13) 



and is given by 



B 3 a 2 - (b + l/2f 
gK) 2(a-6) 2 (2a + 6) 



(4.14) 



for < / < 2(a — b). A graph is presented in 
figure 22. 




Figure 22 The probability density Vg(l) for 
a solid ball when the isometry g is a (pure) 
reflection, eq.( ^4.14p . Here a = 2 and 6 = 1. 
The integrated area is 1. □ 



5 Discussions 



We initially aimed to write out one single ex- 
pression for the probability density Vg(l) for 
screw motions of solid balls, valid for whatever 
values of the four parameters a,b,t , and uo. 
However, we soon found that such expression 
would demand a quite large number of step 
functions to account for all sort of possibilities. 
Since in practice the isometries are dealt with 
one at each time, we found more appropriate 
to present a simple method to have the exact 
Vg(l) for each individual screw motion with 
fixed values of the four parameters. Neverthe- 
less, for those isometric motions of solid balls 
described by three or less free parameters the 
exact expression for Vg(l) for any euclidean 
isometry is short enough and was displayed. 

As promised in the Introduction, we ex- 
hibited the analytic counterpart Vg(l) of the 
computer simulations of pair separations his- 
tograms of the euclidean isometries in cosmic 
crystallography thus far obtained in the liter- 
ature. The graph of Vg(l) in figure 10 corre- 
sponds to the isometries b and c in the Fagun- 
des and Gausmann |J study, or equivalently 
the isometries (3 and S in Gomero J7|. The 
discontinuity in I ~ 0.7 is not observed in the 
two upper figures 1 (Universe E4) of due to 
the strong statistical noise present in these his- 
tograms; nevertheless it is clearly seen in the 
position s = I 2 = 0.5 in the mean histogram 
5b of 0| as well as in the position / ~ 0.7 
in the mean histograms la, lb, 2b, 4a, and 
5 of [H]. Similarly, figure 11 corresponds to 
the isometry a of || and a in and 0; the 
discontinuity in I — 1 has their counterparts 
again in the histogram 5b of [[J and in the 
histograms la, lb, 2a, 4a, and 5 of [§]. 

In contrast with the screw motions (figures 
10 and 11), the pure rotations (figures 13, 15, 
and 16, all with t = 0) do not show disconti- 
nuity of Vg(l) ■ Oppositely to figure 13, where 
a < b , figures 15 and 16 correspond to a > b , 
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so the ball B now has fixed points and the References 



graph of Vg(l) effectively starts from I = 0. 
The strange-looking graph in figure 15 was 
confirmed in a computer simulation; the ir- 
regular behavior near I — 0.7 corresponds to 
the narrow r-interval where the pair-of-discs 
combined intersection eq.(|3.48|) occurs. 



In figure 16 we have chosen values for a, b, 
and u) such that the three types of combined 
intersection (ring, pair of discs, and one disc) 
have equal range in the / scale. En passant, the 
pair-of-discs /—range is now wide, and does 
not originate a bump as did in figure 15. 

Figure 17 corresponds to rotation of the 
solid ball B around a diameter T>; from ( p. 51 ) 
we find that defining l max = 2asina;/2 then 
the graph of l m axPg{l) against I /l max does not 
depend on l max . The points of B along the di- 
ameter T> are fixed under the isometry, so the 
graph again effectively starts from the origin. 

Figure 20 corresponds to a glide reflection 
whose sections x = const in the intersection 
B fl B g are either bicircular discs (for small 
\x\) or circular (for larger The mini- 

mum displacement I occurs for the points of 
B in the intersection with the reflector plane 
x = 0, giving l min = t, the translation. Since 
for all points near the reflector plane we have 
I ~ i(l+2|x| 2 ), then these points are displaced 
almost the same value I ~ t; as a consequence, 
Vg(l) diverges in the vicinity of I = t. Never- 
theless the integrated area is finite, with value 
1. 

Between I = t = 0.25 and I ~ 0.5 we have 
bicircular sections (|4.6| ), while for 0.5 < / < 
lmax ~ 2.0 the sections BnB g nX x are circles. 

The transition from a glide reflection to a 
(pure) reflection is worth describing: if in fig- 
ure 20 we continuously displace the vertical 
line I = t = 0.25 towards 1 = 0, then the 
region of divergence of V®{1) shrinks continu- 
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ously and disappears when / 
giving the graph of figure 22. 



0, eventually 
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